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ABSTRACT 

Wc present optical and near-infrared observations of the dim afterglow of GRB 
020124, obtained between 2 and 68 hours after the gamma-ray burst. The burst 
occurred in a very faint {R > 29.5) Damped Lja Absorber (DLA) at a redshift of 
z = 3.198±0.004. The derived column density of neutral hydrogen is log(iVj^7) = 
21.7 ± 0.2 and the rest-frame reddening is constrained to be E(B — V) < 0.065, 
i.e., Av < 0.20 for standard extinction laws with Ry ~ 3. The resulting dust-to- 
gas ratio is less than 11% of that found in the Milky Way, but consistent with 
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the SMC and high-redshift QSO DLAs, indicating a low metallicity and/or a 
low dust-to-metals ratio in the burst environment. A grey extinction law (large 
Rv), produced through preferential destruction of small dust grains by the GRB, 
could increase the derived Ay and dust-to-gas ratio. The dimness of the afterglow 
is however fully accounted for by the high redshift: If GRB 020124 had been at 
z = lit would have been approximately 1.8 mag brighter — in the range of typical 
bright afterglows. 

Subject headings: cosmology: observations — dust, extinction — galaxies: abun- 
dances — galaxies: ISM — gamma rays: bursts 



1. Introduction 

Spectroscopy of the optical afterglows of cosmological GRBs allows detailed studies of 
the chemical and kinematical properties of gas along the lines of sight. Independently of the 
brightness of the host galaxy, the resulting insight into the GRB environment and its dust and 
gas content can provide clues to the nature of GRB progenitors as well as the properties of 
high-redshift galaxies. In particular, constraints on the column density of neutral hydrogen, 
N(H I), can be obtained through Lja absorption if the burst is sufficiently distant to redshift 
the Lya line into the near-UV/optical domain. Unlike X-ray studies where one infers N(H I) 
only indirectly based on an assumed metallicity, the derived N(H I) will be independent of 
metallicity. 

In this regard, spectroscopy of GRB optical afterglows bears some resemblance to that of 
Damped Lya Absorbers (DLAs, Wolfe et al. (1986)), which are gas-rich absorption systems 
intervening the fines of sight to background sources, such as QSOs. DLAs are chemically 
enriched and are therefore most likely caused by gas in, or close to, galaxies (e.g., Lu et 
al. (1996) and references therein; Pettini et al. (1997a)). In some cases galaxy counterparts 
of DLAs have been detected directly (M0ller & Warren 1993, 1998: Djorgovski et al. 1996; 
Fynbo, M0ller & Warren 1999; M0ller et al. 2002). The important and interesting difference 
between QSO DLAs and GRB absorption systems is that the lines of sight through QSO 
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Programme 165.H-0464) 
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DLAs are absorption cross-section selected whereas GRB absorbers (at the GRB redshift) 
are 'GRB progenitor site' selected. Therefore, GRB absorbers are likely to probe sightlines 
through different parts of DLA galaxies than QSO absorbers do, and hence contribute to a 
more complete picture of the properties of high-redshift galaxies. 

The dust content of DLAs has been constrained mainly by two methods: i) by comparing 

the colors of background QSOs having intervening DLAs with the colors of a control sample 
of QSOs without intervening DLAs (Pei ct al. 1991), and ii) by studying the metal abundance 
ratios, most notably Zn vs. Si and iron group elements, and inferring the dust depletion by 
comparison with interstellar clouds in the Galaxy (Lu et al. 1996; Kulkarni, Fall & Truran 
1997; Pcttini ct al. 1997b; Ledoux, Bergeron & Petitjean 2002). The main conclusion of these 
studies is that DLAs in general are metal poor and that the fraction of the metals bound in 
dust grains is roughly the same as, or somewhat lower than, that found in the Galaxy. 

This paper presents optical/near-IR photometry and optical spectroscopy of the after- 
glow of GRB 020124. GRB 020124 (burst trigger on Jan 24.44531 2002 UT) was the second 
HETE-II burst (Ricker et al. 2002) to be localized to arcsec precision and is the only burst for 
which a host galaxy has not been detected when searched for with HST (Berger et al. 2002; 
Bloom, Kulkarni & Djorgovski 2002). We determine the redshift of the burst to be 3.20, 
which is the fourth highest redshift known for a GRB to date, and the third highest (after 
GRB 000131 at ^ = 4.50 (Andersen et al. 2000) and GRB 030323 at ^ = 3.37 (Vreeswijk et 
al. 2003a)) based on afterglow absorption lines rather than host galaxy emission lines. We 
also show that it has a very high column density, qualifying it as a DLA (the highest value 
of N(H I) of known GRB afterglows and higher than almost all QSO DLAs). We discuss 
the imphcations for the dust-to-gas ratio of the GRB surroundings, including the effects of 
metaUicity and destruction of dust by the GRB itself. 

We assume a cosmology where Hq = 70 km s~^ Mpc~^, Qm = 0.3, and = 0.7. For 
these parameters, a redshift of 3.20 corresponds to a luminosity distance of 27.47 Gpc and a 
distance modulus of 47.19. One arcsecond corresponds to 7.55 proper kpc, and the lookback 
time is 11.5 Gyr. 

2. Observations 

The observing log and photometric results are reported in Table 1 and summarized 
below. 
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2.1. Optical imaging 

The position of the burst was observed at RIKEN with the 0.25 m f/3.4 hyperboloid 
astrograph equipped with unfihered CCD camera AP7p. The field of view was 50' x 50' 
which covered the entire 12' radius HETE-II error circle (Ricker et al. 2002). The observation 
started at Jan 24.530 UT; 126 frames of 20-s exposure were acquired by Jan 24.574 UT. PSF 
photometry was applied to each of the 126 frames. The resultant photometric measurements 
were combined to yield a 3. Oct detection at the position of the afterglow reported by Berger 
et al. (2002). 

R-band images of the optical afterglow were obtained with StanCam at the 2.56 m 
Nordic Optical Telescope (NOT) between Jan 26.035 and Jan 26.076 2002 UT. I-band im- 
ages were obtained at VLT-Melipal (UT3) between Jan 26.311 and 26.318 UT. R-band images 
were obtained with FORSl on VLT-Melipal between Jan 27.278 UT and 27.292 UT. The 
latter images have a pronounced background pattern due to instrument refiections of the 
very bright background (the moon was ~ 75 % ilhiminated). The variation in the pattern 
across the field is approximately 3-4 %. To fiattcn each individual image we used SExtrac- 
tor to create a smoothed background image which was subsequently subtracted. The five 
background-subtracted images were then combined. 

All images were preprocessed using standard tools and calibrated using nearby compar- 
ison stars from the field photometry of Henden (2002). The derived magnitudes are listed 
in Table 1. 



2.2. Nectr-infrcired imaging 

Near-IR observations were obtained with ISAAC on VLT-Antu (UTl). The data consist 
of a total of 30 min in the Ks band and 30 min in the Js band obtained on Jan 26.2 UT and 30 
min in the Ks band obtained on Jan 27.2 UT. The reduction was carried out with the "eclipse" 
software package (Devillard 1997) and the IRAF "Experimental Deep Infrared Mosaicing 
Software" xdimsum. Eclipse was used to remove effects of electrical ghosts from science and 
calibration frames and to construct flat flelds and bad-pixel maps from a series of twilight- 
sky flats. Sky subtraction and combination of the dithered science frames were carried out 
with xdimsum. To account for the rapidly varying background, the sky value in each pixel 
was calculated as the running median of the pixel value in the exposures taken immediately 
before and after a given exposure. To minimize the effects on the running median of bright 
objects in neighboring frames, the subtraction was carried out in two iterative steps: First, 
a cosmic-ray cleaned, background-subtracted combined image was constructed from which a 
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mask of all the objects was created. This mask was then registered to the individual frames 
and taken into account in the following running median background subtraction. Residual 
background signatures were subsequently removed by fitting high-order polynomials to the 
rows and columns of the individual (masked) frames. Finally, the individual background- 
subtracted images were registered and median combined. Excerpts of the resulting images 
are shown in Fig. 1. 

Photometric calibration to J and Ks magnitudes was performed using standard star 
observations of 9106/S301-D and 9149/S860-D (Persson et al. 1998) on the nights of obser- 
vation. The images were astrometrically calibrated using comparison stars from the Guide 
Star Catalogue II (GSC-II)^^ using the WCStools software package by D. Mink^^. The abso- 
lute astrometry has an rms of Of! A3 in the Ks band. The resulting coordinates of the infrared 
afterglow of GRB 020124 are R.A. (J2000.0) = 09'^32"^50.82^ Dec. (J2000.0) = -11°31'11'.'0, 
consistent with the position of the radio and optical afterglow (Berger et al. 2002). 



2.3. Spectroscopy 

Spectroscopic observations were obtained on Jan 26.3 with FORSl on VLT-Melipal 
using the 300V-I-10 grism and the order-separation filter GG375 in long-slit mode with a 
1.0 arcsec wide slit. This configuration provides a wavelength coverage 3650 - 7500 A at 
a resolution of 13 A. The total integration time was 1200 sec, evenly divided between two 
exposures. Standard data reduction procedures were performed using MIDAS. The two 
separate exposures allowed us to reliably reject cosmic-ray events before co-addition. The 
spectrophotometric standard LTT9491 (Hamuy et al. 1994) was observed using the same 
configuration. 



^^The Guide Star Catalogue II is a joint project of the Space Telescope Science Institute and the Osserva- 
torio Astronomico di Torino. Space Telescope Science Institute is operated by the Association of Universities 
for Research in Astronomy, for the National Aeronautics and Space Administration under contract NAS5- 
26555. The participation of the Osscrvatorio Astronomico di Torino is supported by the Italian Council for 
Research in Astronomy. Additional support is provided by European Southern Observatory, Space Tele- 
scope European Coordinating Facility, the International GEMINI project and the European Space Agency 
Astrophysics Division. 

^^http: / /tdc-www.harvard.edu/TDC.html 
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3. Results 

3.1. Spectral energy distribution: spectral index and redshift estimate 

The extracted, flux calibrated spectrum is shown in Fig. 2. Due to the faintness of the 
object spectrum and the resulting dominance of systematic errors from sky subtraction close 
to bright airglow lines, the signal-to-noise fluctuates rapidly in the red part of the spectrum. 
What may look like absorption or emission features redwards of 5500 A should therefore not 
be trusted. The strong absorption line around 5100 A is clearly significant, and several other 
absorption features bluewards of 5500 A are marginally significant at the 3-4 a level. 

One or more of the absorption lines seen in the blue part of Fig. 2 could be due to Lya, 
which would indicate a redshift in excess of 2.3. Our first objective was therefore to construct 
the spectral energy distribution (SED) in the optical to determine its slope, and to look for 
possible large-scale signatures (spectral drops) which might indicate the onset of the Lyman 
Forest and/or the Lyman Valley as expected in an object with redshift significantly in excess 
of 2.3. For this purpose we carefully defined a series of intervals along the spectrum where 
systematic errors from night-sky emission-line subtraction did not pose any problems. The 
total number of counts was determined in each interval of the afterglow spectrum as well as 
of the spectrum of the spectrophotometric standard. Uncertainties were calculated directly 
via propagation of photon statistics. 

The resulting fluxes and corresponding uncertainties (after correction for airmass) in 
each bin are plotted in Fig. 3. In Fig. 3a the minimum-x^ pure power- law fit (/,^ oc v~^\ 
13 — 2.36 ± 0.23; reported uncertainties are la errors throughout this paper) is overplotted, 
but is seen to be a very poor fit to the data points. In particular the data point at 5100 A 
falls many standard deviations below any power-law fit. Such a large drop in that wide 
a bin can effectively only be caused by a damped Lya line at a redshift of about 3.2. At 
this redshift one expects to see not only the Lya line but also a significant drop due to the 
Lyman Forest, and the signature of the onset of the red slope of the Lyman Valley (M0ller 
& Jakobsen 1990). In Fig. 3b we therefore plot the same data points, but here we overlay 
a model power law including the predicted absorption due to the intervening intergalactic 
medium &X, z — 3.2. The power law again represents the minimum-x^ fit and has an index 
oi P — 1.32 ± 0.25. Comparing the of the fits in Fig. 3a and Fig. 3b, we find that even if 
we ignore the presence of the Lya line (i.e. consider only the large-scale SED), the z = 3.2 
model SED represents an acceptable fit while the pure power law is rejected at the 98.6 % 
confidence level (Ax^ = 16.4 for 20 degrees of freedom (dof)). 
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3.2. Redshift and H I column density of the Lyct line 

Following the determination of the underlying power-law index and approximate red- 
shift, we proceeded to fit the Lja line itself. The part of the spectrum useful for line- fitting 
is shown again in Fig. 4, here normalized by division with the fitted power-law. Overplotted 
are model absorption spectra with z — 3.198. The model spectra include the Lyman series, 
Si II, Si III, and O I lines. On the blue side of the Lya line we see the expected hne-blanketing 
by the Lyman Forest. At this resolution, lines in the Lyman Forest are therefore not useful 
for a precise redshift determination, and we shall rely only on the Si IIA1260 and O IA1302 
lines at 5294 A and 5468 A for this purpose. The first line is an unblended Si II A1260 line 
providing a best-fit redshift of 3.200 ± 0.004. The O IA1302 line is blended with a much 
weaker Si IIA1304 line and provides a redshift of 3.196 ± 0.004. The combined redshift of 
these two lines, including wavelength calibration errors, is 2;metai = 3.198 ± 0.004 which we 
shall adopt as the systemic redshift of the absorber. 

Because of the line-blanketing on the blue side of the Lya line, only the red part of 
the line profile could be used to constrain the H I column density of the absorber. Within 
the uncertainty of the redshift we found that acceptable fits could be obtained in the range 
log{NHi) = 21.7 ± 0.2. The estimated la range is marked in Fig. 4 as dotted fines. This 
column density is among the very highest values observed for QSO DLAs (Storrie-Lombardi 
& Wolfe 2000). The absorption system is unlikely to be "intervening" (i.e. unrelated to the 
GRB host) because of the very small probability of such a high-column-density intervening 
absorber. More importantly, a significantly higher redshift of the GRB is ruled out by the 
lack of even stronger Lyman Forest absorption redward of 5200 A (Fig. 3). 

3.3. Characteristics of the optical decay 

In addition to the data presented here, we include the data reported by Berger et al. 
(2002) to constrain the R-band lightcurve. In Fig. 5 we plot all the R-band data points vs. 
the logarithm of time since the GRB. The power-law decay index derived from a weighted 
fit to all the data is a = 1.64 ± 0.03 {f^ oc (t — tcRB)"")- This fit is formally rejected with a 
reduced x^/dof = 26/14. Excluding the late-time HST data we derive an acceptable fit with 
a decay index of a = 1.49 ± 0.04 and x^/dof = 6/12. Fitting only to the data prior to 1 day 
after the burst we derive a decay index of a = 1.45 ± 0.06 with x^/dof = 5/10. These results 
confirm the conclusion of Berger et al. (2002) that the lightcurve is becoming progressively 
steeper (a 'break') (see Fig. 5). 
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3.4. Optical/neELT-infrcired SED 

Assuming a power-law decay index a = 1.49 ± 0.04 we can interpolate the RIJKs-band 
data taken around Jan. 26 UT to a common epoch at Jan 26.20 2002 UT. After dereddening 
the contemporaneous RIJKs-band points for Galactic extinction (Schlegel, Finkbeiner & 
Davis 1998), the colors are R-Ks^ 2.93 ± 0.26 and J - Ks ^ 1.33 ± 0.20. These colors 
are consistent with the ones expected of a GRB afterglow (see the shaded region of the 
color-color diagram in Fig. 2 of Gorosabel et al. (2002a)). 

We have used the intrinsic spectral index obtained from the VLT spectrum (/3 = 1.32 ± 
0.25, see § 3.1) to obtain a mock B-band point by a power-law extrapolation of the R-band 
measurement. In order to infer information on the extinction {Ay, E(B— V)), the BRIJKs- 
band SED was fitted with a functional form f^, oc x lO^^-^^", where is the extinction 

in magnitudes at rest-frame frequency u. Ay was parametrized in terms of Ay using the 
three extinction laws reported by Pei (1992), i.e., for the Small Magellanic Cloud (SMC), 
Large Magellanic Cloud (LMC) and for the Milky Way (MW). For comparison purposes we 
also considered the unextinguished case (pure power law spectrum given by jy oc v~^^. 

As shown in Fig. 6 and Table 2 (third column) the fits of the three extinction laws are 
consistent with the data. The derived extinctions are consistent with the no-extinction case 
(in which case /? = 0.91 ±0.14). Based on these fits we can set an upper limit of E(B — V) < 
0.18. This upper limit however corresponds to an unrealistic value of /9. Constraining the 
spectral index to be /5 > 0.5 and the extinction law to be that of the SMC (which provides 
the best fit to this and most other afterglow SEDs) we find that E(B - V) < 0.065. The /3 
values derived for the three extinguished cases agree with the pure power-law case. 

4. Discussion 

4.1. Jet— wind fireball model 

Within the framework of the afterglow synchrotron model (Sari, Piran & Narayan 1998), 
the decay index, a, and the spectral index, /3, are related through the slope of the electron 
energy distribution, p. Given the low extinction inferred from the SED analysis we find 
that a jet [y < z/g, the synchrotron cooling frequency) expanding into a wind medium 
(n oc r~^) provides the best fireball model for the data. For «i = 1.49 ± 0.04, this model 
requires p = 2.32 ± 0.05, and f3 = 0.66 ± 0.03, consistent with the observational results with 
modest extinction. Expansion into a homogeneous medium (z/ < Uc) is not strictly ruled 
out but requires an unusually high value of p: the predicted values are p — 2.99 ± 0.05 and 



P — 1.00±0.03 (for CKi = 1.49±0.04), consistent with the observational values in the absence 
of extinction. 



4.2. How robust is the redshift determination? 

The spectrum has low resolution and signal-to-noise ratio, but there are five independent 
pieces of evidence for the redshift. First, assuming a pure power-law for the spectrum 
provides a spectral index of 2. 36 ±0.2 3 which is incompatible with /3 = 0.91±0.14 determined 
from the photometry. Second, the strong absorption line at 5100 A cannot be explained by 
anything else than Lya at a redshift around 3.2. Third, assuming z = 3.2 and fitting a model 
with free spectral index but now including the predicted Lyman Forest and Valley absorption, 
one obtains a minimum-^^ fit for a spectral index of 1.32 ± 0.25 which is fully compatible 
with the photometric determination. Fourth, the z = 3.2 fit leads to an improvement in 
over a single power-law alone (for the same degrees of freedom) which rejects the single 
power-law at the 98.6 % level. Fifth, there are two identified metal lines, and none of the 
lines predicted from the explicit fit (e.g. Ly/3) are inconsistent with the observed spectrum 
(Fig. 4). 



4.3. Dim or bright? 

The fact that GRB 020124 was a fairly dim burst (Berger et al. 2002), but fitted by 
a normal fireball model without excessive extinction reinforces the conclusion of Hjorth et 
al. (2002) (GRB 980613), Berger et al. (2002) (GRB 020124), and Fox et al. (2003) (GRB 
021211) that there is a population of dim, unextinguished GRB afterglows that can account 
for at least some of the large fraction of ('dark') bursts for which no optical afterglow is 
detected (Fynbo et al. 2001a). In the case of GRB 020124 this is mainly due to the fairly 
high redshift (unlike GRB 980613 and GRB 021211 which are at redshifts 1.096 (Djorgovski 
et al. 1999) and 1.006 (Vreeswijk et al. 2003b), respectively): If an afterglow is redshifted 
from zi to Z2 it is dimmed by A/x-|-2.5(/3 — a — 1) log ((1 + Z2)/{1 + zi)). For a median GRB 
redshift of 2^1 = 1 and Z2 — 3.198 the relative dimming amounts to 1.82 mag (for a = 1.49 
and P — 0.91). If the afterglow had been 1.82 mag brighter it would have been a typical 
bright afterglow: R ^ 16.7 after two hours and R 20.7 after one day (see e.g. Fig. 3 of 
Gorosabel et al. (2002b) or Fig. 2 of Fox et al. (2003)). 
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4.4. N(H I) vs. E(B-V) 

The Galactic relation between the column density of neutral hydrogen and the reddening 
is N(H I)/E(B-V) = 4.93 ± 0.28 x lO^^ cm^ mag"^ (Diplas & Savage 1994). For the LMC 
and SMC the corresponding ratios are N(H I)/E(B— V) = 2 ± 0.5 x 10^^ cm^ mag^^ and 
N(H I)/E(B-V) = 4.4 ± 0.7 x 10^^ cm^ mag-^ respectively (Koornneef 1982; Bouchet et 
al. 1985). In Table 3 we list this ratio for known GRBs for which the H I column density 
and E(B— V) have been measured or constrained from optical spectroscopy and optical/near- 
infrared photometry (the values given are for an SMC extinction law). In addition to GRB 
020124 these are GRB 000301C (Jensen et al. 2001) and GRB 000926 (Fynbo et al. 2001b,c). 
We exclude GRB 021004 which has several intervening absorbers at different redshifts. It is 
evident that the values for N(H I)/E(B— V) are not larger than expected for galaxies like the 
LMC or SMC. Taken at face value these results indicate that the GRB surroundings are low 
in dust, because of low metallicity (e.g., for a universal dust-to-metals ratio) and/or because 
of a low dust-to-metals ratio (a large fraction of the metals being in the gas phase). This 
conclusion is fully consistent with what is found in QSO DLAs (Pettini et al. 1997b). 

It has been suggested that the intense UV and X-ray flux from a GRB may photoionize 
the surrounding gas and destroy dust grains (Waxman & Draine 2000; Pruchter, Krolik 
& Rhoads 2001; Perna & Lazzati 2002), resulting in changes in the derived N(H I), Ay, 
E(B— V), and their ratios. Galama & Wijers (2001) studied the inferred absorption of soft 
X-rays in BeppoSAX GRB afterglows and compared it to the optical extinction as inferred 
from a fireball model fit to the available data. They claimed evidence for a high N(H I) /Ay 
ratio and concluded that dust destruction may be the cause. The most convincing case 
for dust destruction along these lines was made by Galama et al. (2003) who found a high 
value of N(H I) /Ay for GRB 010222, based on Chandra X-ray Observatory data. At the 
same time they interpreted a strong (i.e., non-grey) far- ultraviolet (restframe wavelength 
~ 1100 A) component in the extinction of GRB 010222 as evidence for dust destruction. 

In a comparison of the properties of DLAs and GRB absorbers, Savaglio, Fall & Fiorc 
(2003) claimed that GRB 990123 {z = 1.60), GRB 000926 {z = 2.04), and GRB 010222 
{z = 1.48) (of which only GRB 000926 has a measured N(H I) (Fynbo et al. 2001b)) have 
higher metal column densities and contain more dust than DLAs at similar redshifts. They 
further suggested that the relatively low reddening (restframe wavelength > 1200 A) observed 
towards these GRBs may be due to a grey extinction law, produced by dust destruction 
(Perna, Lazzati & Fiore 2003). 

Castro et al. (2003) argue that the chromium-to-zinc ratio towards GRB 000926 indi- 
cates that the host is depleted in dust relative to local values to a similar degree as QSO 
DLAs at the same redshift. This conclusion is consistent with what we have suggested above 
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for GRB 020124, namely that the GRB surroundings have a low dust-to-gas ratio, similar 
to QSO DLAs (Pettini et al. 1997b). An additional argument in favor of this interpretation 
is that a significant amount of grey extinction would make the afterglow intrinsically even 
brighter than indicated above (§4.3). For example, assuming a Galactic value of N(H I)/ Ay, 
log N(H l)= 21.7, destruction of 2/3 of the dust producing the extinction in the rest-frame 
V band (Perna, Lazzati & Fiore 2003) and a fiat extinction curve would dim the afterglow 
by 1 mag. Any non-grey component would significantly increase this value. 

We conclude that the high value of the ratio between column density and optical extinc- 
tion first found by Galama & Wijcrs (2001) from X-ray spectroscopy remains when N(H I) is 
estimated from optical spectroscopy in GRB afterglows with a DLA host and the extinction 
is estimated from the observed reddening of the optical/NIR afterglow. The dust destruction 
interpretation originally proposed by Galama & Wijers (2001) to account for the X-ray result 
could conceivably also be apphed in this case. The observed small reddening would be due 
to a modification of the extinction law due to preferential destruction of small grains by the 
prompt UV and X-ray radiation from the GRB, leading to a grey extinction law. However, 
we have found that the straightforward alternative explanation, namely that GRBs occur 
in less chemically enriched environments, similar to those of the SMC and QSO DLAs, is 
fully consistent with the optical observations reported here. Future joint optical and X-ray 
spectroscopy of 1.6 < 2; < 2.5 GRBs combined with accurate multi-wavelength observations 
of the afterglow, in particular in the optical/NIR regime, can be used to distinguish between 
these two possibilities and at the same time fix the metallicity of the burst environment. 
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Fig. 1.— Ks images obtained with ISAAC on 2002 Jan 26.20 UT (left) and 27.20 UT (right) 
of GRB 020124 marked by an arrow. The images are 15" x 15". North is up and east is to 
the left. 
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Fig. 2.— Extracted, flux calibrated FORSl 300V spectrum of GRB 020124 obtained on 2002 
Jan 26.34 UT. The ordinate is fx in arbitrary units. Prom 4000 A to 5500 A the signal-to- 
noise ratio per resolution element (13 A) in the continuum grows from 1.8 to 3.5. Longwards 
of 5500 A the signal-to-noise ratio is completely dominated by sky-subtraction errors and 
therefore varies rapidly from pixel to pixel. In order to obtain the best possible value for the 
spectral slope of the spectrum we have identified "good" sections of the spectrum as detailed 
in § 3.1 and shown in Fig. 3. 
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Fig. 3. — Spectrum of GRB 020124 (Fig. 2) summed in bins free from strong sky lines, (a) A 
pure power-law does not provide an acceptable fit mostly due to the strong absorption line at 
5100 A. This strong absorption feature can only be explained by a damped Lya line close to 
z — 3.2. A power-law fit to the data excluding the bin at 5100 A leads to a spectral slope of 
/3 = 2.36±0.23 incompatible with the optical/near-IR imaging, (b) A power-law fit including 
the effect of the Lyman Forest and Lyman Valley for a source redshift z = 3.198 again 
omitting the bin at 5100 A. The pure power-law fit (a) is rejected at the 98.6 % confidence 
level. The power-law slope of = 1.32 ± 0.25 resulting from the absorbed model fit (b) is 
consistent with the slope derived from the optical/near-IR imaging (/5 = 0.91 ± 0.14). 
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Fig. 4. — FORSl 300V spectrum and model overlay (solid curve) of a 2; = 3.198 damped 
Lya line yielding log(iVjy/) = 21.7 ± 0.2. The estimated la range is plotted as dotted lines. 
The model has been fit to the red wing of the Lya line and the two lines redward of it. These 
metal lines, identified as Si IIA1260 and O IA1302 at 5294 A and 5468 A, were used for the 
redshift determination. For illustrative purposes wc also show the region blucward of the 
Lya line. This region is dominated by Lyman Forest absorption; the data are of very low 
signal-to-noise. The predicted model is also shown but is not fit to the data in this region. 
The data are consistent with the model considering the effect of Lyman Forest absorption. 
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Fig. 5. — The R-band lightcurve based on our data (filled circles) and from Berger et al. 
(2002). The dashed line shows a fit to the early data points only. The dotted line is the 
result of a fit including also the two points from 1 to 3 days after the burst. 
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Fig. 6. — The BRIJKs-band spectral energy distribution of the afterglow on Jan 26.20 2003 
UT. The filled circles represent the RIJKs-band measurements. The filled square is the 
fiducial B-band photometric point obtained by extrapolating the R-band point using the 
intrinsic spectral slope (/3 = 1.32 ± 0.25) derived from the VLT spectrum. The fluxes have 
been corrected for Galactic extinction (Schlegel, Finkbeiner & Davis 1998). The frequencies 
are given in the rest frame for 2; = 3.198. 
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Table 1. Log of observations and photometry of the afterglow of GRB 020124 



Date Filter^Grism Exposure time Telescope^ FWHM^ Brightness*^ 

(Jan 2002 UT) (sec) (arcsec) (mag) 



24.55 unfiltered 126 x 20 RIKEN 2.2 18.5;^j 

26.04 R 6 X 600 NOT 1.14 23.42 ± 0.25 

26.20 Ks 30 X 60 Antu 0.50 20.55 ± 0.13 

26.23 Js 20 X 90 Antu 0.68 21.94 ± 0.16 

26.32 I 2 x 60 + 120 Mehpal 0.74 22.94 ±0.15 

26.34 300V 2 x 600 Mehpal 1.2 

27.20 Ks 30 X 60 Antu 0.56 21.921°;^^ 

27.28 R 5 X 180 Melipal 0.66 24.261°;^ 



*The broad-band optical (R, I) filters used were Bessel filters (Bessel 1990). 

''RIKEN represents a 0.25 m f/3.4 hyperboloid astrograph (Wako, Saitama, Japan) 
equipped with an unfiltered CCD camera; NOT is the 2.56 m Nordic Optical Telescope 
(La Palma, Canary Islands, Spain) equipped with StanCam; Antu and Melipal are the 8.2 m 
Unit Telescopes 1 and 3 on ESO's VLT at Paranal Observatory, Chile, equipped with ISAAC 
and FORSl, respectively. 

^Measured seeing full width at half maximum of point-like objects in the field. 

"^Uncorrected for Galactic extinction. The brightnesses are reported as Rc, Ic, J, Ks 
magnitudes. 



Table 2. SED fitting results 



Extinction law xVdof E(B-V) (3 



MW 1.22 0.045±0.101 0.69±0.50 

LMC 0.56 0.098±0.079 0.29±0.52 

SMC 0.31 0.071±0.051 0.31±0.43 

No extinction 0.78 0.91±0.14 
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Table 3. Column density and reddening for GRB absorbers 



GRB 


z 


log N(H I) 


E(B-V) 
(mag) 


N(H I)/E(B-V) 
(10^^ cm^ mag-i) 


Reference 


000301C 

000926 

020124 


2.040 
2.038 
3.198 


21.2 ±0.5 

21.3 ±0.2 
21.7 ±0.2 


0.031 ±0.014 
0.062 ±0.020 
< 0.065 


ri+113 

09+21 

> 49 


Jensen et al. 2001 
Fynbo et al. 2001b,c 
this paper 



